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A C®  (penalty)  finite  element  is  developed  for  the  equations  governing  the 
heterogeneous  laminated  plate  theory  of  Yang,  Norris  and  Stavsky.  The  YNS 
theory  is  a generalization  of  Mindlin's  theory  for  homogeneous,  isotropic 
plates  to  arbitrarily  laminated  anisotropic  plates  and  includes  shear 
deformation  and  rotary  inertia  effects.  The  present  element  can  also  be 
used  in  the  analysis  of  thin  plates  by  appropriately  specifying  the  penalty 
parameter.  A variety  of  problems  are  solved,  including  those  for  which  (ove  ') 
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A PENALTY  PLATE-BENDING  ELEMENT  FOR  THE  ANALYSIS 
OF  LAMINATED  ANISOTROPIC  COMPOSITE  PLATES 


J.  N.  Reddy 

School  of  Aerospace,  Mechanical  and  Nuclear  Engineering 
The  University  of  Oklahoma,  Norman,  OK  73019 


SUWARY 

A C®  (penalty)  finite  element  Is  developed  for  the  equations  governing 
the  heterogeneous  laminated  plate  theory  of  Yang,  Norris  and  Stavsky. 

The  YNS  theory  Is  a generalization  of  Mindlln's  theory  for  homogeneous. 
Isotropic  plates  to  arbitrarily  laminated  anisotropic  plates  and  Includes 
shear  deformation  and  rotary  Inertia  effects.  The  present  element  can 
also  be  used  In  the  analysis  of  thin  plates  by  appropriately  specifying 
the  penalty  parameter.  A variety  of  problems  are  solved.  Including  those 
for  which  solutions  are  not  available  In  the  literature,  to  show  the 
material  effects  and  the  parametric  effects  of  plate  aspect  ratio,  length- 
to-thlckness  ratio,  lamination  scheme,  number  of  layers  and  lamination 
angle  on  the  deflections,  stresses,  and  vibration  frequencies.  Despite 
Its  simplicity,  the  present  element  gives  very  accurate  results. 
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INTRODUCTION 

Over  the  past  few  years  composites,  especially  fiber-reinforced  laminates, 
have  found  Increasing  application  In  many  engineering  structures.  The 
fiber-reinforced  composites  possess  two  desirable  features:  one  Is  their 
high  stiffness-to-welght  ratio,  and  the  other  Is  their  anisotropic  material 
property  that  can  be  tailored  through  variation  of  the  fiber  orientation 


and  stacking  sequence>-a  feature  which  gives  the  designer  an  added  degree 
of  flexibility. 

Recent  developments  In  the  analysis  of  plates  laminated  of  fiber- 
reinforced  materials  Indicate  that  the  thickness  effect  on  the  behavior  of 
the  plate  Is  more  pronounced  than  In  Isotropic  plates.  The  classical  thin 
plate  theory  assumes  that  normals  to  the  midsurface  before  deformation  re- 
main straight  and  normal  to  the  midsurface  after  deformation,  Implying  that 
transverse  shear  deformation  effects  are  negligible.  As  a result,  the  free 
vibration  frequencies,  for  example,  calculated  using  the  thin  plate  theory 
are  higher  than  those  obtained  by  the  Mindlln  plate  theory^,  which  Includes 
transverse  shear  and  rotary  inertia  effects;  the  deviation  Increases  with 
Increasing  mode  number.  Higher  order  linear  theories  that  include  trans- 
verse shear  effects  have  also  appeared  (see  Relssner^  and  Lo,  Christensen, 
and  Wu^).  Elasticity  solutions  by  Pagano  and  his  associates'*"’  Indicate  the 
Inadequacy  of  the  classical  laminated  plate  theory  (e.g.,  Relssner  and 
Stavsky®,  Dong  et  al.®,  and  Bert  and  Mayberry^®,  In  which  the  classical 
Kirchhoff-Love  kinematic  assumptions  are  adopted  and  effects  of  transverse 
shear  deformations  are  neglected.  The  transverse  shear  deformation  effects 
are  even  more  pronounced,  due  to  the  low  transverse  shear  modulus  relative 
to  the  In-plane  Young's  moduli.  In  the  case  of  filamentary  composite  plates. 
A reliable  prediction  of  the  response  characteristics  of  high  modulus  com- 
posite plates  requires  the  use  of  shear  deformable  theories. 

A number  of  shear  deformable  theories  for  laminates  have  been  proposed 
to  date.  The  first  such  theory  for  laminated  Isotropic  plates  Is  due  to 
Stavsky''.  The  theory  has  been  generalized  to  laminated  anisotropic  plates 
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by  Yang,  Norris,  and  Stavsky^^.  a review  of  various  other  theories,  for 
exanjple,  the  effective  stiffness  theory  of  Sun  and  Whitney^^,  the  higher- 
order  theory  of  Whitney  and  Sun^**,  and  the  three-dimensional  elasticity 
theory  of  Srinivas  and  Rao^s,  can  be  found  in  the  paper  by  Berths.  Other 
approximate  theories  that  have  been  proposed  in  the  literature  include  the 
refined  laminated  plate  theory  of  Mau^^,  the  continuum  theory  of  Sun, 
Achenbach,  and  Herrmann i®,  and  diffusing  continuum  theory  of  Bedford  and 
Stern^®  which  were  primarily  developed  for  use  in  wave-propagation  problems. 
It  has  been  shown,  see  for  example,  the  papers  by  Sun  and  Whitney^®  and 
Srinivas  and  Rao^®,  that  the  Yang-Norris-Stavsky  (YNS)  theory  is  adequate 
for  predicting  the  overall  behavior  such  as  transverse  deflections  and 
natural  frequencies  (first  few  modes)  of  laminated  anisotropic  plates. 

The  first  application  of  the  YNS  theory  is  apparently  due  to  Whitney 
and  Pagano^®,  who  considered  cylindrical  bending  of  antisymmetric  cross- 
ply  and  angle-ply  plate  strips  under  sinusoidal  load  distribution  and  free 
vibration  of  antisymmetric  angle-ply  plate  strips.  Fortier  and  Rosettos^^ 
analyzed  free  vibration  of  thick  rectangular  plates  of  unsymmetric  cross- 
ply  construction  while  Sinha  and  Rath22  considered  both  vibration  and  buck- 
ling for  the  same  type  of  plates.  Recently,  Bert  and  Chen^®  presented, 
using  the  YNS  theory,  a closed-form  solution  for  the  free  vibration  of 
simply  supported  rectangular  plates  of  antisymmetric  angle-ply  laminates. 

While  considerable  effort  has  been  expended  in  the  finite-element 
analysis  of  isotropic  plates,  only  limited  investigations  of  laminated 
anisotropic  plates  can  be  found  in  the  literature.  Pryor  and  Barker^**,  and 
Barker,  Lin  and  Dara^^  used  the  conventional  displacement  finite-element 
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method  to  analyze  thick  laminated  plates.  The  element  has  seven  degrees 
of  freedom  (three  displacements,  two  rotations,  and  two  shear  rotations) 
per  node.  Exploiting  the  symmetries  exhibited  by  anisotropic  plates,  Noor 
and  Mathers^^"^®  studied  the  effects  of  shear  deformation  and  anisotropy 
on  the  accuracy  and  convergence  of  several  shear-flexible  displacement 
finite  element  models  based  on  a form  of  Reissner's  plate  theory.  The 
analysis  was  limited  to  symmetrically  laminated  cross-ply  plates  and  the 
element  used  involved  80  degrees  of  freedom  per  element.  The  conventional 
finite  element,  when  applied  to  relatively  thick  laminated  plates,  either 
has  failed  to  predict  accurately  the  loeal  deformations  and  stresses  of  a 
plate  under  bending  or  is  too  expensive  to  use  due  to  the  large  number  of 
degrees  of  freedom  involved  for  even  relatively  simple  problems.  Mau  and 
Witmer  ^^and  Mau,  Tong,  Pian^°  have  employed  the  so-called  hybrid-stress 
finite-element  method  to  analyze  composite  plates  including  shear  deforma- 
tion. The  hybrid  elements  have  proven  (see  Gallagher^^)  to  have  some  con- 
vergence problems,  and  in  some  cases  they  give  erroneous  results.  Most 
recently.  Panda  and  Natarajan^^  used,  following  Mawenya  and  Davies^^,  the 
quadratic  shell  element  of  Ahmad,  Irons  and  Zienkiewicz^**  with  the  same  nor- 
mal rotation  through  the  thickness  to  claim  improved  accuracy  over  Mawenya 
and  Davieses.  The  'thickness  concept'  mentioned  there  is  essentially  the 
same  as  that  used  in  the  YNS  theory  12.  The  authors  were  primarily  concerned 
with  the  accuracy  of  the  element,  and  no  attempt  was  made  to  solve  new  prob- 
lems for  which  there  do  not  exist  any  closed-form  or  exact  solutions. 

The  present  paper  is  concerned  with  the  development  of  a simple  C® 
element  for  YNS  theory  of  laminated  composite  plates.  The  penalty  function 


placements  and  two  slopes  (i.e.  shear  rotations),  per  node.  The  accuracy 
of  the  element  Is  demonstrated  via  problems  for  which  the  exact  solutions 
and  numerical  results  are  available,  and  results  are  also  presented  for  a 


variety  of  problems  for  which  solutions  are  not  available  in  the  literature. 


LAMINATED  PLATE  THEORY  OF  YAN6-N0RRIS-STAVSKY  (YNS) 

Consider  a plate  of  constant  thickness  h composed  of  a finite  number,  L,  of 
thin  anisotropic  layers  oriented  at  angles  0i,  02,  ...,  0j_.  The  origin  of 
the  coordinate  system  Is  located  within  the  middle  plane  (x-y)  with  the 
z-axis  being  normal  to  the  mid-plane.  The  material  of  each  layer  is  assumed 


to  possess  a plane  of  elastic  symmetry  parallel  to  the  xy-plane.  We  shall 
denote  the  middle  plane  with  q. 

The  YNS  theory  Is  based  on  the  following  assumed  displacement  field: 
u * Ug(x,y,t)  + zi<i^(x,y,t) 

V » v^(x,y,t)  + zi/y(x,y.t)  (1) 

w » w(x,y,t) 

i ; 

where  u,  v,  and  w are  the  displacement  components  in  the  x,  y,  and  z i, 

I . 

directions,  respectively,  t is  the  time,  u and  v are  the  in-plane  |' 

0 0 j 

(stretching)  displacements  of  the  middle  plane,  and  and  ^ are  the  i 

X y 

shear  rotations.  Recalling  the  strain-displacement  equations  of  linear  | 

elasticity,  we  have  j 
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3u  3V  an, 

2.  + 2 — e = 2.  + 2 e =0 

3x  3x  ’ y ay  ^ ay  ’ 
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^xy  ay  ax  ' ay  ax ' • 
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Owing  to  the  existence  of  a plane  of  elastic  symrvetry,  the  constitutive 
relations  for  any  layer  in  the  (x,y)  system  are  given  by 
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where  Q. . are  the  (material)  stiffness  components. 

^ J 

Introducing  the  stress  and  moment  resultants  per  unit  length. 
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we  can  write  (2)  and  (3)  in  terms  of  the  resultants  and  displacements: 
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The  material  componets  and  D^j  are  given  by 
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The  stiffness  coefficients  Q^j  depend  on  the  material  properties  and 


orientation  of  the  m-th  layer.  The  parameters  are  the  shear  correction 


coefficients. 

The  equations  of  motion  associated  with  YNS  theory  are 

3 = 4'., 
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where 
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0^“^  being  the  material  density  of  layer  m,  and  P * P(x,y)  is  the  trans- 
versely distributed  load. 

The  strain  energy  and  the  kinetic  energy,  for  a fixed  time  t,  are 
given  by 

U « Ui  + U2 
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Note  that  the  quantities  in  the  square  brackets  of  U2  are  the  shear  forces 


Q and  Q , respectively. 
X y 


PENALTY  FUNCTION  FORMULATION  OF  THE  EQUATIONS 
The  assumption  of  the  classical  thin-plate  theory  that  the  normals  to 
the  midsurface  before  deformation  remain  straight  and  normal  to  the  mid- 
surface after  deformation  implies  that 


, 3w 

“ • dflO  ® 

3x  ’ ''  '^y  3y 


(13) 


If  we  substitute  for  -ii  and  ^ from  (13)  into  (11),  we  obtain  the  strain 

X y 


energy  U = Uj  associated  with  the  classical  thin-plate  theory.  In  that 
case  U|  involves  the  second-order  derivatives  of  the  transverse  deflection, 
and  the  associated  (conventional)  finite-element  formulation  results  in 
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complicated  elements  (with  many  degrees  of  freedom).  Approaches  that  have 
been  taken  to  relax  the  continuity  requirements  placed  on  the  shape  functions 
in  the  displacement  formulation  of  the  thin-plate  theory  include,  in  addi- 
tion to  the  nonconforming,  hybrid  and  mixed  formulations , the  "discrete 

37  38 

Kirchhoff  hypothesis"  of  Wempner,  Oden  and  Kross  , Fried  , and  the 

39 

"residual  energy  balancing"  and  "reduced  integration"  techniques  of  Fried  , 

. . . ‘♦0  41 

Zienkiewicz,  Taylor  and  Too  , and  Hughes,  Taylor  and  Kanoknukulchai  . The 

present  penalty  function  method  is  a formalization  and  extension  of  these 

ideas  to  the  shear  deformable  theory  of  laminated  composite  plates. 

The  problem  of  finding  the  static  solution  (u,v,w)  to  the  thin  plate 

equations  can  be  viewed  as  one  of  finding  the  critical  points  of  the  total 

potential  energy  = Uj  + V,  where  U;  is  the  strain  energy  given  by  (11) 

in  terms  of  u,  v and  w,  and  "V  is  the  potential  energy  due  to  applied  loads. 

Alternately,  the  problem  can  also  be  viewed  as  one  of  finding  (u,v,ii»  ,4>  ) 

A y 

subject  to  the  constraint  conditions  in  (13).  To  incorporate  the  constraints, 
one  can  use  the  Lagrange  multiplier  method,  or  the  penalty  function  method. 

If  the  Lagrange  multiplier  method  is  used,  we  have 


\ = Ui 


[a  (-^  + ) + X (-^  li'  )]  dxdy 

X 3x  X y 3y  y ■'  ^ 


(14) 


where  A and  A are  the  Lagrange  multipliers.  Comparing  the  Euler 
X y 

equations  of  + V with  those  of  ir  = U + V,  we  see  that  the  Lagrange 

multipliers  are  given  by. 


X = Q = A44  (t^  ■*■’(>)■'■  A45  (-^^  4*  ) 

X ^x  '3x  ^x'  '3y  ^y' 

X = Q = A45  (-1^  +4*  ) Acc  (•—  + 4*  ) 

y ^y  3x  '3y  ^y' 


(15) 
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Thus,  Is  equivalent  to  the  strain  energy  U of  the  shear  deformable 
theory. 

If  the  penalty  function  method  is  used,  the  modified  functional  is 
given  by  ir^  ■ Uj  + Up  + V,  wherein  the  penalty  functional  Up  is  chosen  to 
be 

“p  ■ 7 ’ll  * * ‘i  ’W  * */  * <|7  * ‘xXl?  * 

where  and  the  penalty  parameters.  Clearly,  in  the  limits  ei,  e2-*  “, 

the  constraints  are  satisfied  exactly.  As  opposed  to  the  Lagrange  multiplier 
method  the  constraints  are  satisfied  only  approximately,  and  no  additional 
variables  are  introduced  In  the  penalty  method.  Comparing  the  Euler  equations 
of  the  functional  Hp  with  the  equations  of  the  YNS  theory,  we  see  the 
correspondence, 

^x  ^x  rax  12  'ay  ^y^ 


Qy*»Qj  = ejGj  (|J+  + g2(0+ 


2 ay 


gJ  » Ai,4  , , g^  « ^55 

This  correspondence  implies  that  for  very  large  values  of  g^ , the  equations 
govern  the  thin-plate  theory,  and  for  values  of  g^  given  in  (17),  the 
equations  coincide  with  the  YNS  theory. 

FINITE-ELEMENT  MODELS 

Here  we  present  a (semi discrete)  finite-element  model  based  on  iv  (u,v,\i)  ,w), 

r ^ y 

We  assume,  over  each  element  n^,  the  same  kind  of  interpolation  for  all  of 
the  variables, 

u*  ■ S u?  N?  , V*  ■ S vf  N?  , etc.  (n  ■ nodes  per  element)  (18) 
0^11  0^11 
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where  N®  are  the  element  interpolation  (or  shape)  functions,  and  u®,  and 
V®  are  the  nodal  values  of  u®  and  v®,  respectively.  Substituting  (18) 
into  the  first  variation  of  ir®  (u,v,w,i(/^,i(»y) , and  collecting  the  coefficients 
of  the  variations,  6u^,  6v^,  etc.,  we  obtain 

[M®]{?}  + [K®]{A®}  = {F®}  (19) 

where 


{U®} 

pCs°]^ 

R[S°]^ 

{V®} 

"p[S°] 

"rCs°] 

{W®} 

. [M]  = 

'pCs°] 

(»') 

\ 

symmetric  I[S°] 

(♦J) 

The  elements  (a,6=l ,2, . . . ,5;  i ,j*l  ,2,. . . ,n)  of  the  stiffness  matrix, 
S°j  of  the  mass  matrix  are  given  by 

Kj!  ■ A„  Sjj  » A,e  (S-[5  * * As, 

qj  • * A,e  Sjj  * Aae  * A,e  s]>( 

K|5  ■ 8„  S!|j  ♦ B,e  (SJ5  ♦ Sj()  * S^j 
KiJ  . Bi2  + Bi6  sjj  + 026  S^[j  + Bee 

"ij  • * A,,  sJfj  . Aes 

k;3  ■ B,s  Sfj  * Bes  Sjj  * B,a  * B,e  S^j 
KJI  - B,,  (Sjj  * Sjj)  * Bse  Sjj  * B^a  S^fj 
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S ^1j  *^2  ^1j  *1®2 


'‘ij  • M s”  • «i'i  s{5  . i^l!  ■ 'i'2  5“  ♦ .!  s{° 

»tj  ■ “uSjj  ‘ D,e(Sj}  * Sj{)  ♦ 0„sJ[j  > ej  S»J 
K?]  • 0,2S^  . De.sJ^  ♦ O.eSjj  * * c,e,  S»j 

Kf]  • 0,,(S«J  * Sjp  * OesS^j  ♦ DjjSjj  * .|  sjj 


1(13  ■ 1(23  m 0 ■ 


P N 


^ dxdy  , F|  « F2  • F^f  « F|  - 0 


(21) 


For  free  vibration,  equation  (19)  becomes 


([K®]  - a.2[M®])  {A®}  « {0} 


(22) 


where  w is  the  frequency  of  the  natural  vibration.  For  static  analysis, 
{A}  is  set  to  zero.  The  element  stiffness  matrices  are  assembled  in  the 
usual  manner,  and  boundary  conditions  of  the  problem  are  imposed  before 
solving  for  {a}  or  oj^^. 

In  the  present  study  linear  (n«4)  and  quadratic  (n»8)  elements  of 
the  serendepity  family  are  ussd.  The  element  stiffness  matrices  for  these 
elements  are  of  order  20x20,  and  40x40,  respectively. 
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NUMERICAL  EXAMPLES  AND  DISCUSSION 
The  (penalty)  finite  element  developed  herein  was  employed  in  the 
bending  and  free  vibration  analyses  of  a variety  of  layered  composite 
rectangular  plates.  All  of  the  numerical  results  presented  here  were 
obtained  using  a uniform  mesh  of  2x2  quadratic  (i.e.  8-node  quadrilateral) 
elements  in  the  quarter  plate.  Computations  were  conducted  on  an  IBM  370/ 

158  computer  in  single  precision. 

In  the  following  analyses  two  types  of  boundary  conditions,  simply 
supported  and  clamped,  and  two  types  of  orthotropic  materials  were  used. 

The  coordinate  system  and  boundary  conditions  are  shown  in  Figure  1.  The 
material  properties  used  are  (0^2  * v,,  * '>'13)  • 

Material  I:  E1/E2  “ 25,  “ 0.5,  S23/E2  “ 0.2,  '^32  “ 0.25. 

Material  II:  E1/E2  » 40,  012/^2  * 0.6,  G23/E2  = 0.5,  'Jh  “ 0.25 

A value  of  5/6  was  used  for  the  shear  correction  coefficients  (see  Whitney**^). 


y, 


Figure  1 Coordinate  system,  finite-element  mesh,  and  boundary 
conditions 
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Bending  Analysis 

First  the  effect  of  the  reduced  integration  on  the  bending  deflection 
and  stresses  is  examined  using  a four-layer,  cross-ply  (0°/90°/90®/0®) , 
simply  supported  square  plate  (Material  II)  subjected  to  sinusoidal  loading 
(SSL), 

P * P-  cos  ^ cos  ^ 

0 a b 

The  percentage  error  (between  the  solution  obtained  by  using  2x2  Gauss 

points  and  3x3  Gauss  points)  in  the  center  deflection  and  maximum  normal 

stress  (a  s a ) as  a function  of  the  side-to-thickness  ratio  (a/h)  are 
A y 

shown  in  Figure  2.  The  stresses  were  computed  at  the  Gaussian  points 
using  equation  (3).  Figure  3 shows  the  bending  deflection  versus  the  side- 
to-thickness  ratio  for  the  same  problem  using  2x2  Gauss  rule.  This  result 

43 

is  in  excellent  agreement  with  the  closed- form  solution  of  Whitney  . 

Thus,  the  standard  3x3  Gauss  rule  (for  the  numerical  integration  of  ele- 
ments in  equation  (21))  gives  less  accurate  results,  especially  for  ratios 
a/h  > 10.  Guided  by  this  observation  (also,  see  Zienkiewicz  et  al.‘*o) 
the  remaining  results  were  obtained  using  2x2  Gauss  rule. 

Figure  3 also  shows  the  stresses,  a , o , and  T for  the  four-layer, 

X y xy 

cross-ply  (0®/90°/90°/0°) , simply  supported  square  plate  under  sinusoidal 
loading.  To  see  the  effect  of  loading  and  material  on  the  deflection,  the 
same  problem  was  solved  using  Material  II  and  uniform  loading,  and  Material 
I and  sinusoidal  loading.  Note  that  decreasing  the  ratio  E1/E2  from  40  to 
25  has  the  same  effect  as  using  the  uniform  loading  in  place  of  sinusoidal 
loading.  Bending  deflections  and  stresses  are  presented  in  Table  1 for  a 
(4-ply  (0“/90‘’/90°/0®) , Material  I)  clamped  plate  under  sinusoidal  loading, 
and  simply  supported  plate  under  point  load  at  the  center. 


iaterial  1 1 , SSL , SS) 


w,  SSL,  Material  I 


w,  SSL,  Material  II 


(0V9OV9OVO^  SS) 


(SSL,  Material  II) 


w,  UOL,  Material  II 


a/h  “ 10 


Figure  3 Bending  deflections,  and  stresses  vs.  side  to 
thickness  ratio  for  four-layer,  cross-ply, 
simply  supported  square  plates 


7.  (SSL,  Material  II) 
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To  further  illustrate  the  accuracy  of  the  present  element,  two  problems 
for  which  exact**^  and  finite  element  solutions^^’^^  available,  were 
solved  and  results  are  sunmarized  in  Tables  2 and  3 (for  Material  I).  Table 
2 contains  the  normalized  bending  deflections  and  stresses  for  three-layer, 
cross-ply  (0“/90*/0®),  simply  supported  square  plate  subjected  to  sinusoidal 
loading.  The  outer  layers  are  each  h/4,  and  the  middle  layer  is  h/2  thick 
(i.e.  sandwich  construction).  Table  3 contains  similar  information  for  three - 
layer  (equal  thickness),  cross-ply  (0®/90®/0“),  simply  supported  rectangu- 
lar (b/a  ■ 3)  plate  under  sinusoidal  loading.  Present  solutions  are  com- 
pared with  exact  solutions  of  Pagano^  and  Pagano  and  Hatfield  , and  the  finite- 
element  solutions  of  Panda  and  Hatarajan  "^and  Maweny  and  Davies  . It  is  clear 
that  the  present  solution  is  the  closest  to  the  exact  solution  for  the 
deflection  for  all  ratios  of  a/h.  Since  the  stresses  in  the  present  study 
are  computed  at  the  Gaussian  points,  it  is  not  meaningful  to  compare  for 
relative  accuracy. 

Figure  4 shows  the  normalized  bending  deflections  versus  the  ratio 
a/h  for  the  3-layer,  cross-ply,  simply  supported  square  plate  under  sinus- 
oidal loading.  For  Material  I the  same  problem  was  resolved  with  layer 
orientation  of  0®/91*/0®  (the  middle  layer  is  now  oriented  at  91°)  to  see 
the  effect  of  slight  variation  (introduced,  say,  in  manufacturing)  in  the 
orientation  of  the  layers  on  the  deflection.  Note  that  the  error  in  the 
angle  causes  slight  variation  in  the  deflection  only  at  higher  values  of 
a/h  (i.e.  for  thin  plates). 

Figure  5 shows  plots  of  bending  deflection  versus  the  side -to -thickness 
ratio  for  twy-layer,  cross-ply  (0°/90°)  square  plate  (Material  II)  under 
sinusoidal  and  uniform  loadings,  and  for  four-layer,  symmetric  angle-ply 


a 


I 


I 
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Table  2 Three-layer  (0'’/90®/0®)  simply  supported  square  plate  subjected  to 
sinusoidal  loading  (Material  I,  tj  “ ta  » h/4,  t?  ■ h/2) 


Normalized 

center 

a/h 

Source 

deflection 

v7 

Normal  stresses,  o (top  and  bottom)* 


-Oy(0.0.h/4) 


/a  b h\ 


12.5  Present  FEM 


Pagano  & Hatfield^ 
Present  FEM 
Panda  & Natarajan^^ 
Mawenya  & Davies 
Present  FEM 


Pagano  & Hatfield^ 
50  Present  FEM 

Panda  & Natarajan^^ 

Mawenya  & Davles^^ 

Pagano  & Hatfield^ 
100  Present  FEM 

Panda  & Natarajan^^ 

Mawenya  & Oavles^^ 


Classical  plate  theory 


1.189 

1.1364 

1.114 

1.273 

1.0866 


1.031 

1.0197 

1.016 

1.048 

1.008 

1.0055 

1.003 

1.015 


1.000 


0.543 

0.5118 

0.557 

0.546 

0.5154 


if  ■ Wa(h3/PQa‘*)  , 7 • ah^/p^a^  , a »{4Gi2  + [Ei  + (l+vi2)E2]/(l-vi2V2i)}  it'*/12 
♦Computed  at  the  Gaussian  points  In  the  present  study  . 
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Table  3 Three-layer  (0®/90®/0®)  simply  supported  rectangular  plate 
(b/a  ■ 3)  subjected  to  sinusoidal  loading  (Material  I) 


Normalized 

center 

deflection, 

w 

Normalized 

stress,  o'  (top  and  bottom)* 

a/h 

Source 

+o^(0.0,h/2) 

5;o,(0.0,h/6) 

5 

Present  FEM 

1.695 

0.5984 

0.0691 

0.01789 

6.25 

Present  FEM 

1.267 

0.6006 

0.0540 

0.01338 

Exact:  Pagano^ 

0.919 

0.725 

0.0435 

0.0123 

10 

Present  FEM 

0.802 

0.6031 

0.0364 

0.01017  1 

Panda  & Natarajan^^ 

0.752 

0.653 

0.0367 

0.0105 

Mawenya  & Davies 

1.141 

0.685 

" 

0.0141  ! 

12.5 

Present  FEM 

0.694 

0.6038 

0.0322 

0.00941  i 

1 

Exact:  Pagano^ 

t 0.610 

0.650 

0.0299 

0.0093  j 

1 20 

Present  FEM 

0.578 

0.6045 

0.0276 

0.00858  i 

' 

i 

Panda  & Natarajan^^ 

0.565 

0.654 

0.0287 

0.0091 

i 

Mawenya  & Davies 

0.664 

0.651 

0.0099 

25 

Present  FEM 

0.551 

0.6046 

0.0264 

0.00838  i 

I 

j 

Exact:  Pagano^ 

0.520 

0.628 

0.0259 

0.0084  1 

150 

Present  FEM 

0.515 

0.6044 

0.0251 

0.00812  I 

Panda  & Natarajan^^ 

0.513 

0.654 

0.0264 

0.0087 

Mawenya  & Davies^^ 

0.529 

0.640 

0.0087  ! 

Exact:  Pagano^ 

0.508 

0.624 

0.0253 

0.0083  1 

1 100 

Present  FEM 

0.506 

0.6034 

0.0253 

0.00802  1 

1 

Panda  & Natarajan^’ 

0.505 

0.654 

0.0261 

0.0086 

Mawenya  & Davies^^ 

0.510 

0.638 

-- 

0.0085 

Classical  plate  theory 

0.503 

0.623 

0.0252 

0.0083 

w ■ w 100  E2h3/P^a‘*,  T * ah^/P^a^ 

* Computed  at  the  Gaussian  points  in  the  present  study* 
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(45°/-45®/-45®/45°)  square  plate  (Material  I and  II)  under  sinusoidal 
loading.  It  is  clear  that  the  effect  of  shear  deformation  is  quite  sig- 
nificant in  (cross-ply,  as  well  as  angle-ply)  composites  with  side-to- 
thickness  ratio,  a/h  < 20. 

Bending  deflections  and  stresses  versus  side- to- thickness  ratio  are 

shown  in  Figure  6 for  two-  and  four-layer  angle-ply  (45°/-45°/45‘’//45°) 

square  plates.  For  the  same  plate  thickness,  a two-layer,  angle-ply  plate 

undergoes  larger  deflection  than  the  four-layer,  angle-ply  plate.  The 

stresses  'a  = a and  ? are  shown  only  for  the  four-layer,  angle-ply 
A y xy 

plate  under  uniform  loading.  The  stresses  corresponding  to  the  sinusoidal 
loading  (for  the  same  problem)  are  scaled  with  respect  to  the  stresses 
associated  with  the  uniform  loading.  In  the  case  of  sinusoidal  loading 
both  normal  stress  and  shear  stress  increase  with  decreasing  side  to  thick- 
ness ratio. 

Free  vibration  analysis 

Figure  7 shows  plots  of  the  nondimensionalized  fundamental  frequencies 
versus  side- to- thickness  ratio  for  4-layer,  angle-ply  (symmetric  and  anti- 
symmetric) square  plates.  The  result  obtained  for  simply  supported  (Material 
II)  plate  is  in  good  agreement  with  the  closed- form  solution  of  Bert  and 
Chen^^.  The  present  study  predicts  higher  frequencies,  with  the  deviation 
increasing  with  a/h.  Figure  7 also  shows  the  plot  of  fundamental  frequencies 
for  the  symmetric  angle-ply  (45°/-45°/-45°/45° , Material  I).  Incidentally, 
this  plot  is  in  excellent  agreement  with  that  in  Figure  5 of  Whitney  and 
Pagano^®.  However,  the  figure  caption  there  (i.e.  in  reference  20)  says 
that  the  result  was  obtained  for  four-layer,  antisymmetric  angle-ply 
(45®/-45®/45®A45°) , simply-supported  square  plate  (Material  II).  As  pointed 


i 


SSL , Material  1 1 
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out  by  Bert  and  Chen^^,  and  confirmed  by  the  present  study,  the  plot  shown 
in  Figure  5 of  Reference  20  does  not  correspond  to  the  antisymmetric  angle- 
ply  plate.  To  identify  the  result  with  the  right  problem  the  author  experi- 
mented with  Material  I and  with  clamped  boundary  conditions,  for  which  results 
are  also  shown  in  Figure  7.  Obviously,  none  of  these  come  close  to  that 
presented  by  Whitney  and  Pagano^o.  Thus,  Figure  5 of  Whitney  and  Pagano-® 
corresponds  to  four-layer,  symmetric  angle-ply  (45°/-45'’/-45°/45'’) , simply 
supported  square  plate  with  layers  made  of  Material  I. 

Similar  results  are  presented  in  Figure  8 for  three-  and  four-layer 
cross-ply  simply  supported  plates  (Material  II).  Interestingly,  the  three- 
layer  and  four-layer  cross-ply  square  plate  have  almost  the  same  funda- 
mental frequencies  for  a/h  < 15.  Figure  9 shows  plots  of  nondimensional- 
ized  fundamental  frequencies  versus  the  angle  of  orientation  for  three- 
and  four-layered  square  plates  (Material  II,  and  a/h  » 10).  Finally, 

Figure  10  shows  the  effect  of  the  number  of  layers  on  the  fundamental  fre- 
quency of  layered  angle-ply  (45V-45°/+/-/. . . ) square  and  rectangular  plates 
(Material  II).  For  number  of  layers  greater  than  six,  the  fundamental 
frequency  is  virtually  the  same. 

CONCLUSIONS 

Using  the  penalty  function  concept  of  Courant^^  the  equations 
governing  the  thin-plate  theory,  a shear  deformable  theory  for  layered 
composite  plates  that  resembles  the  YNS  theory’^  is  presented.  A finite 
element  model  based  on  the  penalty/YNS  theory  is  developed  herein  and 
applied  to  the  bending  and  free  vibration  analyses  of  rectangular  composite 
plates  with  various  edge  conditions  and  loadings.  The  numerical  results 
are  compared  with  those  obtained  by  other  finite-element  methods  and  with 


-o~  (e/-e/a/-e,  Material  II) 
(9/-a/a,  Material  I) 


Figure  9 Fundamental  frequencies  vs.  angle  of 
orientation  for  angle-ply,  simply 
supported  square  plates  (a/h  « 10) 


b/a  ■ J 


2 4 6 8 10  12  14  16  18  20 

No.  of  layers 

Figure  10  Fundamental  frequencies  vs.  number  of 
layers  for  four-layer,  angle-ply 
(45®/-45’/45®/-45®) , simply  supported 
plate  (a/h  ■ 10,  Material  II) 


29 


exact  solutions.  The  present  element,  despite  its  simplicity  in  formulation 
and  programming,  gives  the  most  accurate  results. 

Application  of  the  element  to  nonlinear  (in  von  Karman  sense)  and  bimodu- 
lus (i.e.  different  elastic  properties  in  tension  and  compression)  plate  prob- 

44  »45 

lems  was  investigated  recently  by  the  author  . However,  its  application 
to  a nonlinear,  shear  deformable  theory  of  composite  plates  is  still  awaiting. 
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APPENDIX 

extensional,  flexural-extensional , and  flexural  stiffnesses 
(i.j»l,2,6) 

plate  planform  dimensions  in  x,  y directions 

layer  elastic  moduli  in  directions  along  fibers  and  normal 
to  them,  respectively 

force  components  in  the  finite  element  formulation  (i»l  ,2,. . . ,5)  ’ 
layer  in-plane  and  thickness  shear  moduli 
total  thickness  of  plate 

rotary  inertia  coefficient  per  unit  midplane  area  of  lamina 

shear  correction  coefficients  associated  with  the  yz  and  xz 
planes,  respectively  (i»l,2) 

element  stiffness  coefficients  (i ,j“l  ,2, . . . ,80;  a,3»l ,2, . . . ,5) 
total  number  of  layers  in  the  plate 

stress  couple,  and  stress  resultant,  respectively  (i»l,2,6) 

element  mass  coefficients  (i,j»l,2 8) 

element  shape  functions  {i»l ,2, . . . ,8) 
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n 


nodes  per  element 


P 

P 


P 

0 


R 


laminate  normal  inertia  coefficient  per  unit  midplane  area 
transversely  distributed  load 
intensity  of  transversely  distributed  load 
shear  stress  resultants 

plane  stress  reduced  stiffness  coefficients  (i,j»l,2,6) 

laminate  rotary-normal  coupling  inertia  coefficient  per 
unit  midplane  area 

element  matrices  in  FEM  formulation  (i  ,j»l  ,2, . . . ,8;  i,n«0,x,y) 
time 


displacement  components  in  x,  y.  z directions,  respectively 

in-plane  displacements  in  x,  y directions 

nodal  values  of  displacements  u,  v (i»l,2,...,8) 

strain  energies 

potential  energy 

position  coordinates  in  cartesian  system 

■ shear  strains 

column  of  vector  of  generalized  nodal  displacements 
• penalty  parameters  .(i*1.2) 
normal  strains 

' orientation  of  m-th  laminate  (m*l ,2,. . . ,L) 

■ Lagrange  multipliers 

■ total  potential  energy  functionals 

■ normal  stresses 
> shear  stresses 

■ slope  functions 

■ fundamental  frequency  of  free  vibration 
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